Određivanje uopštene kompleksne funkcije
clear variables
syms w C L R
assume(0 < R & 0 < C & 0 < L )
 
zamena = L == C*R^2
zamena = 
 
vrednosti = [C == 1, R == 1]
vrednosti = 
 
syms s
Z1 = R + 1/(s*C)
Z1 = 
Z2 = 1/(1/(s*L) + 1/R)
Z2 = 
syms H(s)
H(s) = subs(Z2/(Z1 + Z2), lhs(zamena), rhs(zamena))
H(s) = 
 
simplify(H(s))
ans = 
Amplitudska i fazna karakteristika
syms Hjw(w)
Hjw(w) = subs(simplify(H(s)), s, 1i*w)
Hjw(w) = 
syms Aw(w)
Aw(w) = simplify(abs(Hjw(w)))
Aw(w) = 
syms phi(w)
phi(w) = expand(angle(Hjw(w)))
phi(w) = 
Crtanje amplitudske i fazne karakteristike
 
figure
yyaxis left
fplot(w, subs(Aw(w), lhs(vrednosti), rhs(vrednosti)), [0,4])
ylabel('A(w)')
yyaxis right
fplot(w, subs(phi(w), lhs(vrednosti), rhs(vrednosti)), [0,4])
ylabel('\phi(w) [rad]')
xlabel('w')
grid on
Određivanje propusnog opsega
Aref = limit(Aw(w), w, Inf)
Aref = 
assume(0<w)
w3dB = solve(Aw(w) == Aref/sqrt(2), w)
w3dB = 
BandPass3dB = [ w3dB, Inf]
BandPass3dB = 
 
wg1 = subs(w3dB(1), lhs(vrednosti), rhs(vrednosti))
wg1 = 
Označavanje propusnog opsega u grafiku
figure
fplot(w, subs(Aw(w), lhs(vrednosti), rhs(vrednosti)), [0,4])
hold on
rectangle('Position',[double(wg1), 0, 4, double(Aref/sqrt(2)) ], ...
'EdgeColor','r',...
'FaceColor',[1 0.9 0.9],...
'LineWidth',1)
hold off
ylabel('A(w)')
xlabel('w')
grid on